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rj ■ ABSTRACT 

The evolution of number density, size and intrinsic colour is determined for a volume- 
i-C ■ limited sample of visually classified early-type galaxies selected from the HST/ACS images 

of the GOODS North and South fields (version 2). The sample comprises 457 galaxies over 
q 320 arcmin 2 with stellar masses above 3 • 10 10 MQin the redshift range 0.4<z<1.2. Our data 

5— i ■ allow a simultaneous study of number density, intrinsic colour distribution and size. We find 

r/} ' that the most massive systems ( >, 3 • 1O 11 M ) do not show any appreciable change in co- 

& . moving number density or size in our data. Furthermore, when including the results from 

2dFGRS, we find that the number density of massive early-type galaxies is consistent with 
no evolution between z=1.2 and 0, i.e. over an epoch spanning more than half of the current 
age of the Universe. Massive galaxies show very homogeneous intrinsic colour distributions, 
' featuring red cores with small scatter. The distribution of half-light radii - when compared 

\f} , to z~0 and z>l samples - is compatible with the predictions of semi-analytic models relat- 

lO ■ ing size evolution to the amount of dissipation during major mergers. However, in a more 

speculative fashion, the observations can also be interpreted as weak or even no evolution in 
comoving number density and size between 0.4<z<1.2, thus pushing major mergers of the 
most massive galaxies towards lower redshifts. 

| Key words: galaxies: evolution — galaxies: formation — galaxies: luminosity function, mass 

function — galaxies: high redshift 



1 INTRODUCTION tight mass-metallicity relation and abundance rati os suggesting a 

quick build-up of the stellar component (see e.g. lRenzini|[2006h . 

During the past decades the field of extragalactic astrophysics has _ , , . ... , . r , r 

r On the other hand, semi-analytic models of galaxy formation pre- 

undergone an impressive development, from simple models that ,. . , , , , . ^ . r .. , r 

° r . , , diet a more extended assembly history (if not star formation) from 

were compared with small, relatively nearby samples to current . n ,, ,. . , , , . , , 

., major mergers. By carefully adjusting these models, it has been 
surveys extending over millions of Mpc at redshifts beyond z~ 1 ., , , ,. .. ... .... ... , 

J ° r J possible to generate realizations that are compatible wi th the ob- 

along with numerical models that can probe cosmological volumes , . „ ... ... , . , — : — : — ~r\ 

c f . ° served stellar populations in these galaxies (e.g. Kavirai et al. 

with the aid of large supercomputers. However, in the same period |„„„ J |_ T — . — — 71~ nn A IT! — — i IU™A 

. , „ , . , , . , , , , 2006; De Lucia et al. 2006; Bower et al. 2006) 

of time, our knowledge of the baryon physics relating the dark and 

luminous matter components has progressed much slower, mainly m this paper we study the redshift evolution of a sam- 

due to the highly non-linear processes that complicate any ab initio p i e of t he most massi v e earl y-type galaxies from the cata- 

approach to this complex problem. logue of iFerreras et ail d2009h . which were visually selected 

The evolution of the most massive galaxies constitutes one of from the HST/ACS i mages of the GOODS North (HDFN) and 

the best constraints one can impose on the modelling of galaxy for- South (CDFS) fields dGiavaiisco et al. 1120041) . Our data set com- 

mation. Within the current paradigm of galaxy growth in a ACDM plements recent wor k exploring the issue of size and stellar 

cosmology, massive galaxies evolve from subsequent mergers of mass evolution (e.g. |Bundv et alTI 120051: iMcIntosh et all 120051; 
smaller structures. The most massive galaxies are early-type in 



morphology and are dominated by old stellar populations, with a 



Franceschini et al. 2006; Fontana et al. 2006; Borch et al. 200i 
Brown et al. |2007l ; lTruiillo et al. Il2007l : fvan Dokkum et al. M20Q: 



The coverage (320 arcmin 2 ), depth (la surface brightness limit per 
pixel of 24.7 AB mag/arcsec 2 in the i band) and high-resolution 
E-mail: ferreras@star.ucl.ac.uk (FWHM~ 0.12 arcsec) of these images allow us to perform a con- 
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Figure 1. Sample of massive spheroidal galaxies extracted from the v2.0 ACS/HST images of the GOODS North and South fields (Ferreras et al. in prepara- 
tion). We show in black the volume-limited sub-sample used in this paper. We further subdivide this sample into three mass bins, between 3 • 10 10 MQand 
10 12 Mq (stellar mass), and in four redshift bins (top axis) centered around z=0.5, 0.7, 0.9 and 1.0 (corresponding to look-back times of 5, 6.3, 7.3, and 7.7 Gyr 
before present, respectively). The dashed line tracks the characteristic stellar mass (M„) from GOODS-MUSIC Montana et al. |20O6l) . The curved solid lines 
are the limiting (%ab = 24) masses for two e xponentally decaying star formation histories with solar metallicity, started at redshift zp=i and with formation 
timescales r = 1 (top) and 8 Gyr (models of iBruzual & Chariot 2003). Our sample is classified according to a best fit template which roughly separates into 
'red' and 'blue' galaxies, represented in this figure as filled and hollow circles, respectively. 



sistent analysis of the redshift evolution of the comoving number 
density, size and intrinsic colour of these galaxies. 



2 THE SAMPLE 

The HSU ACS images of the GOODS North and South fields (v2.0) 
were used to perform a visual classification o f spheroidal galaxies. 
This is a continuation of lFerreras et al. I d2005h - that was restricted 
to the CDFS field. However, notice that our sample does not apply 
the selection based on the Kormendy relation, i.e. the only con- 
straint in this sample is visual c lassification. The ana lysis of the 
complete sample is presented in iFerreras et al. I j2009h . Over the 
320 arcmin 2 field of view of the North and South GOODS/ACS 
fields, the total sample comprises 910 galaxies down to zab = 24 
mag (of which 533/377 are in HDFN/CDFS). The available pho- 
tometric data - both space and ground-based - were combined 
with spectroscopic or photometric redshifts in order to determine 
the stellar mass content. Spectroscopic redshifts are available for 
66% of the galaxies used in this paper. The photometric redshifts 
have an estimated ac curacy of A(z)/(1 + z) ~ 0.002 ± 0.09 
( IFerreras et aLl [2009). Stellar masses are obtained by convolving 
the synthetic populations of IBruzual & Chariol d2003[) with a grid 
of exponentially decaying star formation histories (see appendix B 



of lFerreras et al. 120091 for details). A lChabrieJ J2003I) Initial Mass 
Function is assumed. Even though the intrinsic properties of a stel- 
lar population (i.e. its age and metallicity distribution) cannot be 
accurately constrained with broadband photometry, the stellar mass 
content can be reliably determined to within 0.2 — 0.3 dex provided 
the adopted IMF gives an accurate represen tation of the true initial 
mass function (see e.g. IFerreras et al. 2008). 

The sizes are computed using a non-parametric approach 
that measures the total flux within an ellipse with semimajor axis 
otot < 1.5ap ct ro- The eccentricity of the ellipse is computed 
from the second order moments of the surface brightness distribu- 
tion. The half-light radius is defined as R50 = %/aso x &50, where 
050 and 650 are respectively the semimajor and semiminor axes of 
the ellipse that engulfs 50% of the total flux. Those values need to 
be corrected for the loss of flux caused by the use of an aperture (see 
e.g. lGrahametal.l2005h . We used a synthetic catalogue of galaxies 
with Sersic profiles and the same noise and sampling properties as 
the original GOODS/ACS images to build fitting functions for the 
corrections in flux and size. The corrections depend mostly on R50 
and, to second order, on the Sersic index. Most of this correction is 
related to the ratio between the size of the object and the size of the 
Point Spread Function of the observations. The dependence with 
Sersic index (or in general surface brightness slope) is milder and 
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for th is correction the concentration (as defined in Bershadv et al. 
2000) was used as a proxy. 

We compared ou r photometry with the GOODS-MUSIC data 
dGrazian et al. 1120061) in the CDFS. Our sample has 351 galaxies 
in common with that catalogue, and the difference between our 
total+corrected i-band magnitudes and the total magnitudes from 
GOODS-MUSIC is Ai = i OUIS - Jmusic = -0.17 ± 0.16 mag. 
This discrepancy is mostly due to our corrections of the total flux. 
A bootstrap method using synthetic images show that our cor- 
rections are accurate with respect to the true total flux to within 
0. 05 mag, and to with in 9% in half-light radius (see appendix A 
of lFerreras et al. Il2009h . Our estimates of size were also compared 
with t he GALFIT-based parametric approach of lHaussler et al. I 
d2007h on the GEMS survey. Out of 133 galaxies in common, 
the median of the difference defined as (Rgg rs -Rg, EMS )/R^ rs is 
—0.01 ± 0.16 (the error bar is defined as the semi-interquartile 
range). 

We focus here on a volume-limited sample comprising early- 
type galaxies with stellar mass M s > 3 x 10 10 M Q . This sample 
is binned according to fixed steps in comoving volume (a standard 
ACDM cosmology with Q m =0.3 and h — 0.7 is used through- 
out). T he complete sample of 910 galaxies from iFerreras et al. I 
is shown in figure [T] Solid (open) circles represent early- 
type galaxies whose colours are compatible with an older (younger) 
stellar population. This simple age criterion is based on a compar- 
ison of the observed optical and NIR colours with the predictions 
from a set of templates with exponentially decaying star forma- 
tion histories, all beginning at redshift zf = 3, with solar metallic- 
ity. The "old" populat ion is compatib le with formation timescales 
t ,$ 1 Gyr (see IFerreras et al.1l2009l for details). The black dots 
in the figure correspond to the sample of 457 galaxies used in this 
paper. 

We further subdivide this sample into three mass bins, starting 
at log(M s /M ) = 10.5 with a width A log(M 3 /M Q ) = 0.5 dex. 
For comparison, the characteristic stellar mass from the mass func- 
tion of the GOODS -MUSIC sample is shown as a dashed line 
dFontana et al~ll2006l) . althou gh we warn that the GOODS-MUSIC 
masses are calculated using a Sa lpeted dl955ri IMF, which will give 
a systematic 0.25 dex overestimate in logM s with respect to our 
choice of IMF. Our sample is safely away from the limit imposed 
by the cut in apparent magnitude (ias 24). The curved solid 
lines give that limit for two extreme star formation histories, cor- 
responding to the "old" and "young" populations as defined above. 
Notice that within our sample of massive early-type galaxies there 
are NO galaxies whose colours are compatible with young stellar 
populations (i.e. open circles). 



3 THE EVOLUTION OF MASSIVE GALAXIES 

The redshift evolution of the comoving number density is shown 
in figure [2] (black dots). The (la) error bars include both Poisson 
noise as well as the effect of a 0.3 dex uncertainty in the stellar mass 
estimates. These uncertainties are computed using a Monte Carlo 
run of 10 ,000 realizations. T he figure includ es data from GO ODS- 
MUSIC dFontana et aLll2006h, C OMBQ17 dBell et al. Il2004r) and 
Pal/DEEP2 dConselice et al. 1120071) . At z=0 we sho w an estimate 
from the segregated 2dFGRS luminosity function dCroton et al. I 
l2005h . We take their Schechter fits for early-type galaxies within 
an environment with a mean density defined by a contrast - mea- 
sured inside radius 8/i _1 Mpc - in the range 6% = —0.43 ■ ■ • +0.32 
(black open circles). In order to illustrate possible systematic ef- 



fects in 2dFGRS, we also include the result for their full volume 
sample as grey open circles. The 2dFGRS data are originally given 
as luminosity functions in the rest-frame bj. We took a range of 
stellar populations typical of early-type galaxies in order to trans- 
late those luminosities into stellar masses. The error bars shown for 
the 2dFGRS data represent the uncertainty caused by this transla- 
tion from light into mass over a wide range of stellar populations 
(with typical M/L(6j) in the range 7 ■ ■ • 12 M / L Q ). The black 
solid lines show semi-a nalytic model (SAM) predictions from 
iKhochfar & Silkl d2006al) . Their SAM follows the merging history 
of dark matter halos generated by the Extended Press-Schechter 
formalism down to a mass resolution of Af m i n = 5 x 10 9 Mq, 
and follow s the baryonic physics with in these halos using recipes 
laid out in Khochfar & Burkertl d2005l and references therein). The 
grey dashe d lines are the pr edicti ons from the Millennium simu- 
lation |Pe Lucia et al. I (e.g. l2006h . This model is extracted from 
their web-based database and is not segregated with respect to 
galaxy morphology. This explains the excess number density in the 
low-mass bin (bottom panel). In the two higher mass bins most of 
the galaxies have an early-type morphology. The predictions of the 
Millennium simulation are in agreement with the middle bin - i.e. 
masses between 10 11 and 3 ■ 1O 11 M0. However, for the most mas- 
sive bin, the sharp decrease in density with redshift of the models 
is in remarkable disagre ement with the observations. In contrast, 
IKhochfar & Silk (2006a) predict a nearly constant density at the 
highest mass bin out to z< 1. 

The main reason for this discrepancy is that AGN feedback in 
the Millennium simulation prohibits the growth of massive galaxies 
by gas cooling and subsequent star formation in order to reproduce 
the right colour-bimodality and the luminosity function at z=0. As 
shown in Khoch far & Silkl d2008l) the existence of a characteristic 
mass scale for the shut-off of star formation will lead to dry merg- 
ing being the main mechanism for the growth of massive galax- 
ies. In that respect the evolution of the number density of massive 
galaxies in the Millennium simulation i s mainly driven by merge rs. 
The difference between that model and lKhochfar & Silkl d2006ah is 
probably due to the different merger rates in their models. The Mil- 
lennium simulation pre dicts a lower major merger rate compared to 
IKhochfar & Silkl ( 2006a ) almost by a factor 10 (Hopkins et al., in 
preparation). 

Figure [3] shows the redshift evolution of the half-light radius. 
Our methodology follows a non-parametric approach avoiding the 
degeneracies intrinsic to profile fitting. Nevertheless, we com- 
pared our size estimates with those using a parametric approach 
like GALFIT (Hauss ler et al. |[2007t) and there is g ood agreement 
(see j 2}. Our data (black dots) are compared with iTruiillo et al.l 
2007, grey triangl es) and with a z~0 measurement from the SDSS 



Shen et al. ]Hjo| taking their early-type sample). The error bars 



give the RMS scatter of the size distribution with in each mass and 
redshift bin. The lines correspond to the models of Khochfar & Silk 
(2006b). These models associate size evolution to the amount of 
dissipation encountered during major mergers along the merging 
history of an early-type galaxy. The points at high redshift (z>1.2) 
correspond to individual measurements from the literature (see cap- 
tion for details). In all the comparisons shown in this paper with 
work from the literature, we have checked that the initial mass 
functions used are similar, so that stellar masses are compa red con- 
sistently. All results quoted either use a IChabrierl (2003J) IMF or 
functions very close to it in terms of the total mass expected per 



1 jhttp : / /www .mpa-garching.mpg. de /millennium! 
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Figure 2. Evolution of the comoving number density of massive early-type galaxies. Our sample (filled dots) is shown for the three mass bins as labelled in the 
three panels. The points and error bars show the average and standard deviation for each redshift bin, adding in quadrat ure Poisso n noise and the uncertainty 
in the stel lar mass estimate. Oth er points from the literature are included, nam ely GOODS-MUSIC (x signs; Fon tana et al. 1200^) . Palomar/DEEP2 (vertical 
segments; Consel ice et al. |2007|) and COMB017 (+ signs; lBorch et al. l|2006j). The points at z=0 are estimates fr om 2dFGRS restrict ed to (spectroscopically 
classified) early-type galaxies (black/grey open circles corresponding to their mean density sample/full volume; ICroton et al. 2005). Notice the density of 
the most massive galax ies (top) does not evolve over the whole redshift range, spanning a look-back ti me of about 8 Gyr. Th e black solid lines are model 
predictions from I Khochfar & Silk 2006a), the grey dashed lines are from the Millennium simulation (e.g. lDe Lucia et al. 120061 ; no morphology segregation) 
and the black dashed lines are fits to a power law: n oc (1 + z) a , with a given in each panel. 



luminosity unit, which - for early-type systems - mainly reduces 
to the shape of the low- mass end of the IMF. Other func t ions u sed 
in the quoted data w ere Kroup a et al. ] Jl993l) ; iKroupal d200ll) or 
iBaldrv & Glazebrookl fe003l) . Only for the GOODS-MUSIC data 
dFontana et aLl l2006 ) the Salpeter IMF (1955) was used, which 
will always give a systematic overestimate of ~ 0.25 dex in stel- 
lar mass with respect to the previous choices given its (unphys- 
ical) extrapolation of the same power law down t o the low stel- 
lar mass cutoff (see e.g. iBruzual & Charlotll2003l) . A single-law 
Salpeter IMF is an unlikely choice for the stellar populations in 
early-type g alaxies as shown by comparisons of photometry with 
kinematics JCappellari et al. Il2006t) or with gravitational lensing 
dFerreras et al. 1120081) . 



Similarly to the density evolution, we also apply a simple 
power law fit only to our data points: R e oc (1 + z) 13 . The solid 
lines give those best fits, and the power law index is given in each 
panel. Taking into account all data points betweeen z=0 and z~2.5 
one sees a clear trend of decreasing size with redshift for all three 
mass bins. However our data suggest milder size evolution for the 
most massive early-type galaxies between z=1.2 and z=0.4, corre- 
sponding to a 4 Gyr interval of cosmic time. 

The depth and high spatial resolution of the ACS images also 
allow us to probe in detail the intrinsic colour distribution of the 
galaxies (i.e. the colour dist ribution within each g alaxy). We follow 
the approach described in iFerreras et al. I {2005) which, in a nut- 
shell, registers the images in the two bands considered for a given 
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Figure 3. Redshift evolution of the half-light radius. Our sample is shown as filled black circle s. The data from Trui illo et al. is shown as grey triangles 

and the z-^0.1 sizes of SDSS early-type galaxies are shown as open circles jShen et al. 120031) . At higher redshif t s we include individual m easurement s from 
Buitrago et al. I i2008i grey '+' signs) and from other recen t work in the literature ( 'x' signs: IZirm et al. Il2007t iToft et al. lEobl : | van Dokkum et al. 1 120081 : 
Cimatti et al. |200S1) . The solid lines are the predictions from lKhochfar & Silkl fe006bl) . and the dashed lines are fits to a power law: R e oc (1 + z)P, with the 
power law exponent given in each panel. 



colour, degrades them by the Point Spread Function of the other 
passband, and perfoms an optimal Voronoi tessellation in order to 
achieve a S/N per bin around 10 while preserving spatial resolution. 
The final binned data is used to fit a linear relation between colour 
and \og(R/R e ) from which we determine the slope and the scatter 
about the best fit (using a biweight estimator). Figure [4] shows the 
observer-frame V— i colour gradient (bottom) and scatter (top) as a 
function of stellar mass (left) and half-light radius (right). The black 
dots correspond to binned data in stellar mass, showing the average 
and RMS value within each bin. Notice the significant trend with 
increasing stellar mass towards redder cores (i.e. more negative 
colour gradients) and small scatter. The colour gradient is in most 
cases nearly flat, and only for the lowest mass bin do we find sig- 
nificantly large gradients. For comparison, we also show as small 
grey d ots a continuation of the original sample from lFerreras et ahl 
(2009) towards lower stellar masses. Blue cores (positive colour 



gradients) dominate in spheroidal galaxies below 10 10 Mq. The 
homogeneous intrinsic colour distribution thereby suggests no sig- 
nificant star formation and a fast rearranging process of the stel- 
lar populations if mergers take place during the observed redshift 
range. Notice this sample only targets objects visually classified as 
early-type galaxies. The early phases of major merging are there- 
fore excluded from our sample. Nevertheless, the number density 
at the massive end (upper panel of figure[2]l does not change signifi- 
cantly between z=0 and z~l, already suggesting that major merger 
events must be rare over those redshifts. 



4 DISCUSSION AND CONCLUSIONS 

Using a volume-limited sample of massive spheroidal galaxies 
from the v2.0 ACS/HST images of the GOODS North and South 
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Figure 4. The radial colour gradient (A(V— i) / A log R; bottom) and RMS scatter (top) is shown as a function of stellar mass (Jeff) and projected physical 
half-light radius (right). Negative (positive) gradients imply galaxies with red (blue) cores. The three mass bins used in this paper are illustrated by the vertical 
lines and the hatched regions in the left-hand panels. The average and RMS values within each bin is shown by the black dots and error bars. Hollow dots in the 
left-hand panels give the values for two bins outside of our adopted mass range, for comparison. The figure shows the trend of a reduced scatter - concentrating 
the distribution towards a slightly red core - as mass increases. 



fields we have consistently estimated the number density, size and 
intrinsic colour distribution over the redshift range 0.4<z<1.2. 
In combination with other samples we find a significant differ- 
ence in the redshift evolution according to stellar mass, in agree- 
ment with recent work based on other samples or different se- 



lection criteria (see e.a.lBundv et al. I 2005; Mcintosh et al. 


2005; 


Franceschini et al. 2006; Fontana et al. 2006; Borch et al. 


2006 


Brown et al. 20ol;lTruiillo et al. 2007: van Dokkum et al. 1 


2008) 



The most massive galaxies - which impose the most stringent con- 
straints on models of galaxy formation - keep a constant comoving 
number density between z~ 1 and (i.e. over half of the current age 
of the Universe) but present a significant size evolution, roughly a 
factor 2 increase between z=l and 0. Note, however that within 
our sample, there is no significant size evolution over the redshift 
range z=0.4- • -1.2. It is by extending the analysis to higher red- 
shifts that the size evolut i on shows up at the m ost massive bin (e.g. 
van Dokku m et al. 1200 8; Buitra go^et al. I200H) . When velocity dis- 
persion is added to the analysis, a significant difference is found in 
the er-R e distribution between z=0 and z=l, su ggesting an impor- 
tant c hange in the dynamics of these galaxies dvan der Wei et al.l 
l2008h . 

Some of the semian alytic models of massive galaxy evolution 
dKhochfar & Silkl2006al lbl) are in good agreement with these obser- 
vations. These models follow the standard paradigm of early-type 
galaxy growth through major mergers, with the ansatz that size evo- 
lution is related to the amount of dissipation during major mergers. 
The decreasing evolution in the comoving number density at high 
masses is explained within the models by a balance between the 
'sink' (loss due to mergers of massive galaxies generating more 



massive galaxies) and 'source' terms (gain from mergers at lower 
mass) over the redshifts considered. One could argue that the sink 
terms would generate a population of extremely massive galaxies 
(above a few 1O 12 M0), possibly the central galaxies within mas- 
sive groups or clusters. However, this population - with predicted 
comoving number densities below 10~ 6 Mpc~ 3 - are very hard to 
study with current surveys. Furthermore, environment effects in 
these systems will compli cate the analysis of size evolution (e.g. 
iKhochfar & Ostrikerll2008h . 

It is important to note that the lack of evolution in the num- 
ber density relates to the bright end of the luminosity function. 
iFaber et al. I d2007» found a significant change in the number den- 
sity of red galaxies with redshift. However, they also emphasize that 
this change does not refer to the most luminous galaxies. If we in- 
clude all mass bins in our sample, we do find a significant decrease 
in the number density with redshift, as the lower mass bins - which 
contribute the most in numbers - do have a rather steep decrease 
in density (see figure [2j. This difference suggests that the (vari- 
ous) mechanisms playing a role in the transition from blue cloud to 
red sequence must be strongly dependent on the stellar mass of the 
galaxies involved. 

In a more speculative fashion, our data are also suggestive 
of weak or even no evolution in the number density of the most 
massive early-type galaxies over a redshift range 0.4<z<1.2. This 
would imply a negligible role of major mergers at the most massive 
end for z>0.4, thereb y pushing this stage of galaxy formation to- 
wards lower redshifts ( Khochfar & Silkl2 008). Another speculative 
scenario for the evolution of massive spheroidal galaxies would in- 
volve negligible major mergers at these redshifts and a significant 
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amount of minor mergers which will 'puff up' the galaxy. Minor 
mergers are considered the cause of re cent star formation o bserved 
in NUV studies of early-type galaxies ( iKavirai etal.l2007h . Larger 
surveys of Luminous Red Galaxies are needed to confirm or dis- 
prove this important issue. 



